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Abstract. The classical problem of thermal explosion is modified so that the chemically 

active gas is not at rest but is flowing in a long cylindrical pipe. Up to a certain section 

the heat-conducting walls of the pipe are held at low temperature so that the reaction rate 

is small and there is no heat release; at that section the ambient temperature is increased 

and an exothermic reaction begins. The question is whether a slow reaction regime will 

>- ' be established or a thermal explosion will occur. 
IT) ! 

The mathematical formulation of the problem is presented. It is shown that when the 
■ pipe radius is larger than a critical value, the solution of the new problem exists only 
up to a certain distance along the axis. The critical radius is determined by conditions 
in a problem with a uniform axial temperature. The loss of existence is interpreted as 
a thermal explosion; the critical distance is the safe reactor's length. Both laminar and 
developed turbulent flow regimes are considered. In a computational experiment the loss 
of the existence appears as a divergence of a numerical procedure; numerical calculations 
reveal asymptotic scaling laws with simple powers for the critical distance. 
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1. Introduction 

The classical problem of a thermal explosion in a long cylindrical vessel containing a 
reacting gas at rest was first formulated and solved by D. A. Frank-Kamenetsky in 1939 
(see ( 1 )). The basic result is that there is no spontaneous thermal explosion if the radius 
of the vessel r is less than a critical value r* determined by the gas properties and the 
properties of the exothermal chemical reaction. If ro > r* , a quiet evolution of the chemical 
reaction is impossible and a thermal explosion occurs. 

In the present paper the following question is asked: What will happen if the long 
cylindrical vessel is in fact a pipe where the chemically active gas is flowing? At a certain 
section of the pipe (see Figure 1) we suddenly increase the temperature of the walls and 
the rate of the exothermic reaction increases. How will the critical conditions for thermal 
explosion be affected by the gas flow? In the present paper we provide a mathematical 
formulation of the problem. Our first result is rather obvious: if the radius of the pipe is 
less than critical at the increased wall temperature, nothing will happen; the reactor will 
be safe for arbitrary length. The second result seems to be deeper: If the radius of the 
pipe is larger than critical, there exists a distance along the pipe up to which the solution 
of the mathematical problem exists, and after which it ceases to exist. This critical length 
depends on the flow conditions and determines the "safe" length of the reactor, up to 
which the thermal explosion does not happen. 

We want to mention here the remarkable paper ( 2 ) where the influence of a flow on the 
critical conditions for thermal explosion was first investigated. The flow under considera- 
tion in ( 2 ) was different-free convection- and therefore the mathematical problem was also 
different from ours. 



2. Problem formulation: Laminar flow 

We make the same basic assumptions concerning the exothermic chemical reaction as 
in the classical theory of thermal explosion (see ( 1 )): Arrhenius dependence of the reaction 
rate on the temperature, large activation energy, and neglect of the longitudinal heat 
transfer in comparison with the contribution of the advection. Then the equation of the 
energy balance takes the form: 

(1) cv(r)d z T = ^d r (rkd r T) + Qa{T)e~ E/RT . 

Here T is the absolute temperature of the gas, c is the heat capacity, z is the longitudinal 
coordinate reckoned from the section where the temperature was elevated, v(r) is the lon- 
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gitudinal fluid velocity, r is the current radius, k is the heat conductivity corresponding 
to the flow regime, Q is the thermal effect of the reaction, a(T) is the pre-exponential 
factor in the Arrhenius reaction rate as a function of the temperature- a slow function of 
the temperature, E is the activation energy, and R is the universal gas constant. Apply- 
ing the Frank-Kamenetsky high activation energy approximation (see ( 1 )), and going to 
dimensionless variables, we reduce equation (1) to the form 

(2) V(p)d c u = -d p p (j^d p <i?J + AV . 

Here p = r/ro is the dimensionless radius, k the temperature diffusivity appropriate for the 
flow regime under consideration, kq the molecular temperature diffusivity, V(p) = v(r)/v, 
v = G/nr 2 , is the average velocity, G the flow discharge rate, Q = z/(G/7ruo)Pr the 
dimensionless longitudinal coordinate, vq the molecular kinematic viscosity, Pr = uq/kq 
the Prandtl number, a characteristic gas constant. Furthermore, A = tq/£, where I = 
(e^ kqTqc/ Q(j)a(Toj) 2 , <p = E/RT , u = E(T — T )/RTq, T being the wall temperature 
at z > 0. For a gas temperature uniformly distributed along the z-axis in a laminar flow 
d z T = 0, k = Ko and we return to the classical problem formulation (see ( 1 )): 

In the case of laminar motion it is easy to show that V(p) = 2(1 — p 2 ), therefore equation 
(2) assumes the form 

(4) (1 - p 2 )d c u = -d p pd p u + AV 

(we redefined C : C = z / (2G / irvo) Pr) . The boundary conditions for the solution u(p, £) 
take the form 

(5) «(P,0) = 0, m(1,C) = 0, lim(p^(p,C)) = 0. 

The first condition (5) reflects the assumption that the temperature of the flowing gas in 
the whole inlet section z = is equal to the temperature of the wall for z > 0. This is not 
essential and can be replaced by an arbitrary inlet condition u(p, 0) = ito(p)- The second 
and third conditions are the same as in the classical problem formulation. It is clear that 
dc_u is positive, therefore for A < A cr = the solution of the problem does exist for 

< ( < 00. For A > A cr the solution to the problem (4)~(5) does exist only at < ( < Co, 
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where Co is a certain finite positive constant which depends on A. The critical dimensional 
length is 

vd 2vro 

(6) -^z = r Q Re Pr Co, Re = — = • 

i/ v 

In section 4 the numerical solution of the problem (4)-(5) will be presented. 



3. Problem formulation: Developed turbulent flow 

In this case we use the power law proposed earlier (see ( 3 )), according to which the 
longitudinal velocity distribution takes the form 



a 

v*y 



3 . \/3 + 5a 



(7) v ( r ) = v ,c(a)l^-) , a=— — , C(a) = 



v Q J 1 2\nRe 1 v ' 2a 

Here y = r$ — r is the distance from the wall, Re = 2vro / vq is the Reynolds number, and 
i>* is the "friction" velocity. 

Equation (2) in this case is reduced to the form 

(8) (1 - p) a d c u = ^d pP 2 (l - pf- a d p u + AV . 

The first two boundary conditions (5) preserve their form in the turbulent case, however 
the third condition — the condition of zero influx at the axis — takes the form 

(9) lim[p 2 (l-p) 1 -^]=0. 

To obtain (8), (9) an expression for the turbulent temperature diffusivity k is needed. This 
quantity is determined by a relation k = A(Re, Pr)vT, where vt is the turbulent kinematic 
viscosity and A(Re, Pr) is a dimensionless factor which we do not necessarily assume to 
be equal to unity. By definition vt = v^p/d y v, so that 

„., A(Re,Pr)Pr ( v.dV'" ,_ a 

The quantity (v^d/v) 1-01 is a function of a, or, which is the same, of the Reynolds number 
Re (see ( 3 )). Substituting this expression into (2), and redefining ( and £ by 

z ■ 2 2a ■ A(Re, Pr) fv*d^~ 2a 

l-a 



^ ^ r aC 2 (a) \ u 
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e ( PK T c \ PrA(Re,Pr) fv*d 



Q(pa(To) J 2 1 -*aC(a) \ u 
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we obtain the basic equation (8). Condition (9) is obtained from (10) and the from the 
condition of zero flux at the axis p = 0. 

The critical dimensional length zq is therefore determined by the relation 



(12) z = r aC 2 (a) [ — 



,\ 2a t 

v*d\ Co 



v J 2 a A(Re,Pr) 

To demonstrate the difference between condition (9) and the corresponding condition 
for the laminar case, consider the equation 

(13) I^ ( i_ p )i-«^ + AV = 0, 

p dp dp 

which, in the subcritical case, describes the temperature distribution at large distances from 
the inlet section, where the longitudinal temperature variation is negligible. If condition 
(9) holds, u(0) is finite, so that integrating (13) we obtain 

(14) Al-^-^-A'jf^ (D^-T^- 

Thus the integral curves of (13) that satisfy condition (9) form a one-parameter family 
having u(0) as parameter. This family possesses an envelope, and the value of the param- 
eter A which corresponds to the intersection of the envelope with the p-axis at p = 1 is the 
critical value. 

In the following section the value of A cr will be determined numerically, in the case of 
turbulent flow, as a function of the Reynolds number, and for values A > A cr the critical 
dimensionless length Co wm also be determined. 



4. Numerical solutions 

It is convenient to use the coordinates p, £ = C/Co m the numerical solution of the 
problems (4)-(5) and (8), (5), (9) so that the domain of integration is fixed: < p < 1, 
< £ < 1; in these coordinates Co appears as a parameter in the differential equations. The 
resulting equations are discretized and yield a system of algebraic equations for u at grid 
points. This system is solved by Newton's method. If Newton's method converges for a 
particular pair A, Co we decide that the pair is inside the existence region in the A, Co-plane. 
If the method does not converge, the pair is outside. 

Newton's method is iterative and an initial condition for iteration must be chosen. If 

this condition is close to an actual solution then the method is rapidly convergent. If the 
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initial condition is not close to the actual solution then the method can fail to converge 
even when the solution does exist. To pick the initial conditions, we first assume that 
for any A, a solution can be found if the value of parameter Co is sufficiently small. We 
pick A, compute a solution for some Co sufficiently small so that the solver converges, and 
solve again starting from the initial conditions given by the preceding solution. When Co 
becomes so large that the solver fails to converge, we say that we have crossed the existence 
boundary. In the laminar case, whenever < A < v^2 5 our method identifies convergent 
solutions for all values of Co tested. When A > \/2, our method identifies convergent 
solutions for all values of Co less than a critical value Co (A). 

Figure 2 shows the graph of Co (A) for the laminar case. The salient feature is the 
divergence of Co (A) at the point A = V2: (o(\/2) = oo. Similar calculations were also 
performed for the turbulent case. Figure 3 shows the graph of A^ r as a function of the 
Reynolds number. Figure 4 shows the graph of the critical values Co(<M- A detailed 
description of the numerical work will be provided elsewhere. 

5. Conclusion 

The formulation of a new mathematical problem for chemically active gas flow in a 
pipe is presented. It is shown that for a supercritical regime, where the radius of the pipe 
is larger than a critical value, the solution of this problem exists only for an interval of 
bounded length along the pipe axis. The loss of the existence of the solution is naturally 
interpreted as a thermal explosion . The dimensional critical length Co is determined both 
for laminar and for developed turbulent flows. 

In the case of laminar flow the results of the numerical computation of Co (A) (A is the 
parameter proportional to the radius of the pipe) are presented in Figure 2. The asymptotic 
branch of the curve Co (A) a t large A corresponds to the scaling law Co (A) = 1.77A -11 ' 4 . 
The results of the numerical computation of A cr as a function of the Reynolds number for 
developed turbulent flow are presented in Figure 3, and the function Co (A) for developed 
turbulent flow at various values of the Reynolds number are presented in Figure 4. An 
interesting property of this graph is the collapse of the curves corresponding to different 
values of the Reynolds number in a very large (3 orders of magnitude) range of Re starting 
from rather moderate values of A. The asymptotic branch of the curve Co (A) a t large A 
corresponds to the scaling law (o(\) = 0.78A -9 / 4 . The origin of the observed scaling laws 
remains unknown to us. 

6 



References 

Zeldovich, Ya.B., Barenblatt, G.I., Librovich, V.B., and Makhviladze, G.M. (1985). 
The Mathematical Theory of Combustion and Explosions (Consultants Bureau, New York, 
London) . 

( 2 ) Kagan, L., Berestycki, H., Joulin, C, and Sivashinsky, G. (1997), Combust. Theory 
Modelling 1, 97. 

( 3 ) Barenblatt, G.I., Chorin, A.J., and Prostokishin, V.M. (1997), Appl. Mech. Rev. 
50, 413. 



7 



Figure Captions 

Figure 1. Chemically active gas flow in a pipe (schematic). 

Figure 2. The dimensionless critical length Co as a function of the basic parameter A 
proportional to the pipe radius, in the case of laminar flow. The critical length tends to 
infinity as A approaches the critical value A cr = a/2 from above, and exhibits the scaling 
law Co ~ A -11 / 4 at large supercritical values of A. 

Figure 3. The critical value of the basic parameter A cr for turbulent flow as a function of 
the Reynolds number Re. 

Figure 4. The dimensionless critical length Co as a function of the basic parameter A in 
the case of developed turbulent flow. Note the collapse of the curves starting from rather 
moderate values of A in a large range of Reynolds number, and the scaling law Co ~ A -9 / 4 
at large supercritical values of A. 
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